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SUMMARY , , r 

Horizon sensors are essential elements in the orientation systems 
of many space vehicles and satellites. The design of horizon sensors, 
as well as an evaluation of their capabilities and limitations, requires 
a knowledge of the spectral radiance emitted by the earth. This report 
investigates the discontinuity between the surface, the atmosphere, and 
outer space for various latitudinal, seasonal, climatic, and meteorologi- 
cal conditions. Calculations are carried out for the following five 
spectral intervals: the water vapor band from 6.33 to 6.85 microns, the 
ozone band from 8.9 to 10.1 microns, the atmospheric window from 
10.75 to 11.75 microns, the carbon dioxide band from 14 to 16 microns, 
and finally the rotational water vapor band from 21 to 125 microns. The 
superior properties of the carbon dioxide band and the rotational water 
vapor region for horizon sensor purposes are shown. 
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THE INFRARED HORIZON OF THE PLANET EARTH* 


by 

R. A. Hanel, W. R. Bandeen, and 
B. J. Conrath 

Goddard Space Flight Center 


INTRODUCTION 

A knowledge of the direction toward the center of the earth is required for the proper orientation 
of many scientific satellites and manned space capsules. Moreover, the required accuracy in deter- 
mining the direction of the local vertical is constantly increasing with the growing sophistication and 
complexity of experiments. 

Normally, a horizon sensor system determines the center of the apparent disk of the earth, which 
is very close to the center of gravity. In order to make the system independent of sunlight, it has 
become customary to restrict the spectral response of the horizon sensor to the infrared part of the 
spectrum or, more specifically, to wavelengths longer than 2 microns. For the sake of simplicity 
and engineering convenience the spectral response has seldom been extended much beyond 15 microns 
(e.g., Reference 1). This utilizes about one- half the energy of a 270°K blackbody, which is very im- 
portant for normally energy limited systems. 

This approach is entirely justified as long as one is satisfied with an uncertainty in the horizon 
on the order of 50 km (about 1 degree for a 600 km orbit) and as long as the system is able to accom- 
modate the changes in signal strength caused by varying meteorological conditions (Reference 2). 
Only recently have horizon sensors that are more selective spectrally been considered. 1 However, 
for the design of advanced stabilization systems, a better understanding of the infrared structure of 
the horizon is required. 

This report investigates the discontinuity between the earth, the atmosphere, and outer space. 
The apparent horizon was computed for several wavelength regions and under a variety of atmospheric 
conditions. Latitudinal, seasonal, and climatic effects were considered, as well as the effect of clouds 
for several altitudes. Calculations are given in detail to benefit the reader who is interested in the 
problem of horizon sensing but who may not be familiar with the problems in the radiative transfer 
of energy. 


•Also published in: jour. Atm. Sci. (March 1963). 
tPrivate communication from D. Goetze, 1962. 
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While writing this paper, the authors discovered that a similar investigation, although for some- 
what different spectral intervals and model atmospheres, had been carried out recently by K. Y. 
Kondratiev and K. E. Yakushevskaya at the University of Leningrad, U.S.S.R. (Reference 3). For ease 
in comparing the two papers, the present results are plotted herein (Figures 10 through 14) in a for- 
mat similar to that used by Kondratiev and Yakushevskaya. 

OUTLINE OF PROBLEM 

In the spectral range of infrared emission, the sharpness of the horizon depends on the vertical 
temperature and pressure profiles and on the distribution of absorbing (and therefore emitting) matter 
in the atmosphere. Water vapor, carbon dioxide, and ozone have strong absorption bands in this 
range, which must be considered. In addition, droplets and ice crystals in clouds, and perhaps dust 
and aerosols, must also be included. 

An attempt has been made to compute the fraction of the radiance N(w/m l 2 -ster) of the earth and 
the atmosphere within particular wavelength regions where only one of the atmospheric constituents 
dominates absorption. Less important constituents have been considered qualitatively only, with the 
exception of the 8.9 to 10.1 micron region, where emission by ozone and the water vapor continuum 
both are taken into account. The wavelength intervals considered are listed in Table 1. 

Table 1 


Wavelength Intervals Considered. 


Wavelength (microns) 

Major Absorber 

Secondary Absorbers 

6.33 - 6.85 

h 2 o 

None 

8.9 - 10.1 

° 3 

H 2 0 continuum 

10.75 - 11.75 

H 2 0 continuum 

Dust 

14 - 16 

co 2 

o 3 ,h 2 o 

21 - 125 

h 2 o 

None 


The intervals between the quoted wavelength regions, such as between 12 to 14 microns and 16 to 
21 microns, are characterized by overlapping bands and semitransparent gaps. The radiance of these 
intervals is somewhere between that of the strong absorption regions of the gases discussed and that 
of the window chosen for this investigation. The ranges less than 6 microns and greater than 125 
microns are of minor interest, since the fraction of thermal energy available in these spectral ranges 
is diminishingly small. 

The fraction of the radiance emitted by the earth and the atmosphere, which can be detected by a 
horizon sensor with a spectral response , is given by 

l r ® rco (9r ( u* ) 

N <* " “ ~ P\ b a < t > du* du* dA. • (1) 

J\=0 J\j~0 
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(All symbols are defined in Appendix A.) For a discussion of radiative transfer equations see, for 
example, Reference 4. In this calculation, the spectral response of a possible detection system has 
been taken as unity within the quoted wavelength ranges and as zero outside. For this specific form 
of < p K we redefine = N Aa in Equation 1. The integration of N Aa was done graphically, with the aid of 
radiation charts. The problem can be divided into the following tasks: 

1. The formulation of an expression for the reduced optical path u*, for a curved atmosphere as 
a function of height, and for various gas pressures, temperatures, and concentrations; 

2. The selection of atmospheres of extreme but realistic properties to show the limits that can 
reasonably be expected in the apparent horizon; 

3. Numerical computation of u* for each gas and atmosphere selected; 

4. Construction of a radiation chart for each spectral range under consideration; and 

5. The plotting and graphical evaluation of the radiation charts for all spectral ranges and 
atmospheres. 

The results of the integrations are summarized in diagrams, and the problem is discussed in the 
sequence previously listed. 

THE REDUCED OPTICAL PATH 

The reduced optical path du* of an absorber has been taken from Elsasser (Reference 5), and its 
relation to the true optical path du is expressed by 

P /T 0 \ 1/2 

du* ” p“ \pf] du . (2) 

The factors P/P 0 and (T 0 /T) 1/2 take into account the pressure and temperature dependence of line width. 
In these computations the above expression for the reduced optical path u* was used for all gases 
except for ozone when the pressure exceeded 10 cm Hg, after which a constant factor Pj/P 0 replaced 
P/P 0 . As pointed out by Elsasser, the results obtained by his replacement are in better agreement 
with the experimental results on ozone absorption given by Walshaw (Reference 6) than are those re- 
taining the P/P 0 factor. 

For all gases at high altitudes, where the absorber concentration is low, most of the lines that 
form an absorption band have not yet reached total absorption in the center of the lines and absorp- 
tion is linearly proportional to the amount of gas and independent of pressure. Therefore, the pres- 
sure and temperature correction applied at these heights introduces a small error. Also, the effect 
of Doppler broadening on lines was neglected but should be considered for heights greater than about 
50 km (Reference 7; see also Reference 8). It is believed, however, that these effects are relatively 
unimportant compared to the uncertainties in the temperature profiles and distributions of water vapor 
and ozone in the higher atmosphere. 
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Figure 1— Geometry for the calculations. 


The computations were performed for a 
concentrically stratified atmosphere (see Fig- 
ure 1)— the distributions of absorbing gases, 
clouds, temperatures, and pressures then 
being functions only of height. 

Since du = q ds, the reduced optical 
path u*, corrected for pressure and temper- 
ature, is given by 


u* = 



f t 0 y /2 

T(s)J q<s)ds 




q(r) dr . 

( r 2 _ ro 2),/2 


(3) 


The integration along the path s has to be performed from the top of the atmosphere to a level r but 
not farther than to the point where the particular ray comes closest to the surface (e.g., for r 0 <R, to 
the surface and, for r 0 > R, to point A in Figure 1). 

The refraction of rays in the density gradient has been neglected. Although the refraction of in- 
dividual rays in the density gradient may not be entirely negligible, in the problem of horizon sensing 
this effect would to the first approximation cancel because of symmetry about the local vertical. For 
this reason, refraction has been neglected in this investigation. For r 0 = 0 we obtain the value of u* 
required for the calculation of the radiance leaving the atmosphere in the direction of the radius vec- 
tor. Equation 3 can be used also to investigate the limb-darkening effect as a function of the zenith 
angle 0 9 since r 0 = r sin 6^. The radiation chart used easily permits the second symmetrical half of 
the atmosphere (to the right of point A in Figure 1) to be taken into account. The radiation chart also 
deals in a simple manner with clouds and the surface of the earth, and therefore these parameters do 
not have to be considered here in the computation of u* . 

The pressure P(h) chosen for a particular model atmosphere was computed separately by numeri- 
cally integrating the hydrostatic equation 


P(h) = P 0 exp- R^ru’) dh ' (4) 

The gas constant for dry air R d has been taken as 2.8704 x 10 6 cm 2 /sec 2 -°K. (The effect of atmos- 
pheric water vapor on the mean molecular weight of air has been neglected; altitude corrections for 
g(h), although small, have been made.) The function q(r) in Equation 3 gives the distribution of the 
absorbing gas. For C0 2 a homogeneous distribution of 0.00031 fractional parts by volume was 


4 



assumed. Consequently, du along the radius vector dr becomes 


du ( cm NTP ) 


dr 


31 


-dr ( km ) 


(5) 


But, since p air P/R ci T, 


du (cm NTP) 


P T n 


31 


dr 


(km) 


and 


u*co f cm NTP) 



( km ) 


Similar considerations for water vapor and ozone yield 


( 6 ) 



For 0 3 , where P/P 0 1 0.132, a constant pressure ratio P l /P 0 = 0.132 was used. The water vapor 
mixing ratio q H 0 is expressed in grams per kg, and q G is given in cm NTP per km. (The distinction 
between specific humidity and mixing ratio of water vapor was neglected in these calculations.) How- 
ever, before the numerical calculations of u* could be carried out, a number of model atmospheres 
had to be selected; this is discussed next. 


THE SELECTION OF MODEL ATMOSPHERES 

The model atmospheres were designed to illustrate the range of radiation characteristics that 
might reasonably be expected in the various spectral ranges considered. Accordingly, emphasis has 
been placed for the most part on the extremes likely to be found in various atmospheric parameters. 
However, in all cases an effort has been made to keep the models plausible and internally consistent. 

Five different temperature profiles were used (Figure 2). These are intended to cover the lati- 
tude variations and, in the case of the arctic atmosphere, the seasonal variations. The profiles are: 

1. Tropical Atmosphere: The tropical profile given in the Handbook of Geophysics (Reference 9) 
was used up to approximately 22 km. Above this level the low-latitude standard profile 
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HEIGHT ABOVE t*tA N SEA LEVEL (* 



TEMPERATURE ( *' K ) 

Figure 2— Model atmosphere temperature profiles. 


proposed by Stroud and Nordberg (Reference 
10), as determined by sounding rocket meas- 
urements on the island of Guam, was used. 

2. Desert Atmosphere: A linearized average of 
two desert soundings was used up to 30 km, 
with the Reference 10 low-latitude standard 
being incorporated above 30 km. This profile 
is typical of low- latitude desert regions. 

3. ARDC Standard Atmosphere (1959): This pro- 
file (Reference 9) is typical of the mid-latitude 
regions. 

4. Arctic Summer Atmsophere: The artic sum- 
mer profile given in Reference 9 was used up 
to 30 km. Above 30 km the high-latitude sum- 
mer profile given in Reference 10 was adopted. 
This combination was chosen to give the 
warmest arctic stratosphere. 

5. Arctic Winter Atmosphere: The arctic winter 
profile as given in Reference 9 was used up to 
24 km. Above 24 km, the high- latitude winter 
profile given in Reference 10 was chosen. This 
combination was adopted to show the influence 
of a cold arctic stratosphere. 


The distributions of the optically active gases were chosen, as were the temperature profiles, to 
illustrate the range of conditions that might be expected. The proportion of carbon dioxide in the at- 
mosphere was taken as constant with respect to altitude and latitude, and as equal to 0.031 percent 
by volume. The water vapor mixing ratios adopted, shown in Figure 3, are described as follows: 


1. Tropical Model: Up to 8 km, mixing ratios based on a linearization of a Trinidad, West Indies, 
sounding (1200 UT, 20 July 1961) were used. This sounding was chosen as representative of 
a very moist tropical troposphere. From 8 to 20 km a constant relative humidity was as- 
sumed. From 20 to 31 km a constant frost point was used, and above 31 km a constant mix- 
ing ratio of 0.11 gm/kg. This combination results in a total of about 5.07 cm of precipitable 
water. Temperature profile no. 1 was incorporated in this model. 

2. Desert Model: A linearized average of two Aoulef, Algeria, soundings (1200 UT, 16 August 
1960; and 1200 UT, 25 August 1960) was used up to 17.5 km. From 17.5 km to 31 km a con- 
stant frost point was assumed, and above 31 km the mixing ratio was held constant at 0.06 
gm/kg. This results in a total of 0.9 cm precipitable water. 

3. Arctic Summer Model: The water vapor mixing ratios were taken from the altitude and lati- 
tude distributions of water vapor given by Manabe and Moller (Reference 11). 
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4. Arctic Winter Model: The mixing ratios for this atmosphere were taken from Reference 11. 

5. ARDC Standard Atmosphere (1959): A representative mixing ratio profile was adopted, re- 
sulting in a total of 1.16 cm precipitable water. 

The four Ozone Models chosen, shown in Figure 4, are described below: 

1. Tropic Model: Temperature profile no. 1 was employed along with the ozone distribution for 
5°N latitude used in Reference 11. 

2. Mid- Latitude Model (Arosa normal): An ozone distribution based on a 3-year mean (1956- 
1958) of measurements made at Arosa, Switzerland, 47 N (Reference 12) was chosen along 
with temperature profile no. 3. 

3. Disturbed Mid-Latitude Model (Arosa high): This model was constructed as an illustration of 
the rapid changes that might be encountered at any given point as the result of the passage of 


frontal systems. The ozone distribu- 
tion used is that given in Reference 12 
based on measurements taken after the 



Figure 3 — Model atmosphere water vapor mixing ratios. 
Above about 15 km, the ARDC and arctic winter mixing 
ratios are coincident with those of the arctic summer. 



Figure 4 — Model atmosphere ozone distribution. 
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passage of a pronounced cold front. As in the Arosa normal case, temperature profile no. 3 
was used. 

4. Arctic Winter Model: Temperature profile no. 5 was used in this model along with the Manabe 
and Moller (Reference 11) distribution of ozone for January at 85 L N. 


NUMERICAL COMPUTATION OF u* 

A large number of values of u* had to be calculated for each case to allow curves of u* versus / 
to be plotted on the radiation diagrams. For these calculations Equations 6, 7, and 8 were solved for 
carbon dioxide, water vapor, and ozone, respectively, on an IBM 7090 Computer. The integration over 
r was carried out by a straight summation using intervals of 0.05 km, initiated at h max = 70 km. 


CONSTRUCTION OF RADIATION CHARTS 

A radiation diagram was constructed for each of the five spectral regions chosen. The various 
types (Moller, Elsasser, Yamamoto, etc.) are similar in their basic concept and can be transformed 
into each other by conformal mapping. Moller f s (Reference 4) version of the diagram and Elasser’s 
(Reference 5) generalized absorption coefficients were adopted. For the computations, each region 
was divided into a number of spectral intervals for which a constant absorption coefficient was as- 
sumed. The values adopted are summarized in Table 2. 


Table 2 


Absorption Coefficients Adopted for the Calculations. 


Major Absorber 

Wavelength K 
(microns) 

Band Absorption, 
r = r (Lu*) 
(Log L) 

Absorption by Continuum, 
T = exp (-ku*) 

(Log k) 

H.O 

6.33 - 6.5 

+ 2 

_ 


6.5 - 6.67 

0.4 + 2 

- 


6.67 - 6.85 

+ 2 

- 

o 3 

8.9 - 9.35 

0.8 - 2 



9.35 - 9.9 

0.2 

_ 


9.9 - 10.1 

0.6 - 1 

- 

Window 




(H 2 0 continuum) 

10.75 - 11.75 

0.7 - 5 

0,85 - 2 

co 2 

14 - 14.5 

0.6 - 1 

_ 


14.5 - 15.5 

0.3 

_ 


15.5 - 16 

0.6 - 1 

- 

h 2 o 

21 - 25 

+ 1 



25 - 30 

0.6 + 1 

_ 


30 - 40 

0.2 + 2 

_ 


40 - 125 

0.6 + 2 

- 
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The radiation diagrams for the carbon dioxide band, the window, the ozone band, and the 6.3- 
micron water vapor band are comparatively simple to construct, since the wavelength regions are so 
small that the distribution of blackbody energy, relative to the spectral intervals into which each re- 
gion has been divided (Table 2), remains approximately constant over the temperature range of inter- 
est. Therefore, lines of constant temperature are horizontal. 

The above simplification, however, does not hold for the broad wavelength region being considered 
within the rotational band of water vapor. For the computation of this radiation diagram, the integral 

B a (T)^d\ % LV' T, A^ ■ ( 9 ) 

Jk i ' 

which is necessary in the construction of the ordinate, was approximated by a summation over the 
same four wavelength intervals shown in Table 2. Here B, A indicates that the actual fraction of the 
Planck function bounded by the somewhat broad limits of each of the four wavelength intervals was 
used, rather than the product of the spectral blackbody radiant emittance at the central wavelength and 
the wavelength interval, which was used in the numerical calculations for all the other regions. 

In the window region, the absorption is mainly determined by the absorption coefficient k. The 
numerical value used in our computation is a factor of 10 higher than the one given by Elsasser (Refer- 
ence 5), which apparently should be used with 10 u*. In all computations, the generalized absorption 
coefficient was considered to be independent of temperature. (The scales in all diagrams were calcu- 
lated with the aid of the radiation slide rule from the Admiralty Research Laboratory, Teddington, 
England.) 


EVALUATION OF THE RADIATION CHARTS 

The atmospheres discussed previously were plotted on the radiation charts, and the integrations 
were performed graphically with a planimeter. In the plotting procedure, the part of the atmosphere 
beyond the closest approach to the earth (to the right of point A in Figure 1) was taken into account by 
imaging the temperature profile in the radiation diagram at the u* value corresponding to r = r Q (Fig- 
ure 5). In cases where the atmosphere was transparent enough to make the surfaces contribute to 
the outgoing radiation, the surface was assumed to radiate as a blackbody at the temperature indicated 
by the lowest point of the particular temperature profile. This effect occurs, for example, in the win- 
dow case for the arctic winter atmosphere (Figure 6). 

In addition to the clear- sky states of the model atmospheres, high cloud overcasts were introduced 
into all except the desert atmosphere. The clouds, too, were assumed to radiate as blackbodies at the 
ambient temperatures associated with the heights of cloud tops. To illustrate extreme cases, the 
cloud- top heights have been considered to be at the tropopause— and therefore at the lowest tempera- 
ture— for the tropical atmosphere (16.8 km), as shown in Figure 7, and for the arctic summer atmos- 
phere (9.0 km). For the standard atmosphere, the cioud-top height has been placed in the minimum 
temperature isothermal region at 12.0 km, a representative level considering the variation in the 


9 



ARCTIC SUMMER, 
ho “ 30 km 


ARCTIC WINTER, h 0 = 30 km 


u* (cm NTP ) 


0 . 

3 0. 

! 

! 


9J 


4 6 8 2 5 

10 1 10 2 


Figure 5—Radiation chart for the 15 micron C0 2 band. Radiance values correspond to the 
area under each curve. The right-hand portion of each curve is an image of the plot to the 
left of point A. The lower arctic winter radiance is caused by lower temperatures and lower 
air densities in the upper atmosphere. 


ARDC, ho = 0 km 



ARCTIC WINTER, h 0 = 0 


{precip. cm 


Figure 6— Radiation chart for the atmospheric window. When the atmosphere is not opaque, 
two solutions exist for h Q = 0, viz., when the ray is terminated at the surface and when the 
ray extends past the rim of the earth. 
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Figure 7— Radiation chart for the 6.3 micron H 2 0 band. The effect of clouds acting as blackbody 
radiators at various heights in the tropical atmosphere Is illustrated. 


height of the tropopause over middle latitudes. Considering the often ill-defined tropopause at high 
latitudes in the winter, the cloud- top height has been placed at 8.5 km for the arctic winter atmosphere. 
This level is near the mean winter tropopause height at high latitudes according to Kochanski (Ref- 
erence 13). 

The radiation chart of the rotational water vapor band (Figure 8) shows samples of the desert at- 
mosphere. The decrease of N (area under each curve) with rising values of h 0 is caused by a rapid 
reduction in the absorbing gas. This rapid decrease is partially compensated by a simultaneous in- 
crease in radiance per unit At as a result of the positive, vertical temperature gradient in the strato- 
sphere and the shape of the isothermal lines in the radiation chart. The point A, for example, rises 
20° from the curve corresponding to h Q = 25 to the one corresponding to 35 km. The high transparency 
of the zone above 25 km shifts the weight of the emitted radiation towards the range that possesses 
the highest absorption coefficient, that is, towards the long wavelength end of the 21 to 125 micron in- 
terval (see Table 2). There the Planck function is less sensitive to small temperature changes, which 
is the physical reason for the rise of the isothermal lines at the left side of the radiation chart. 

Special attention must be given to the second important absorber within each wavelength region. 
Only in the range between 8.9 and 10.1 microns was it necessary to consider two gaseous absorbers 
quantitatively, namely, ozone and the water vapor continuum. This was done, as shown in Figure 9, 
by plotting the calculated transmissivity curve t t • r H ^ 0 and integrating over it. The absorption 
coefficients for the window (10.75 to 11.75 microns) were used to compute t H j0 . 
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2 4 6 8 2468 2 4 68 2 468 

ur 6 io - 5 10 4 io - 3 io - 2 io - 1 io° 


u* (precip. cm HpO) 

Figure 8— Radiation chart for the rotational H 2 0 band. Curves are plotted on the desert 
atmosphere for various values of h 0 . 



Figure 9— Radiation chart for the 9.6 micron 0 3 band. Three curves are illustrated: (1) where 
the transmission is due to 0 3 only; (2) where the first curve is terminated at the equivalent 
blackbody temperature of the window (considering the entire atmosphere); and (3) where the 
transmission is taken as the product of the transmissions due to both 0 3 and H 2 0. Note that 
radiance values calculated from the latter two curves are approximately equal. 
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The influence of water vapor on the 14 to 16 micron carbon dioxide emission band can be neglected 
because of the opaqueness of C0 2 ^.> co 0.02j at a level where water vapor is still transparent 

(; H 0 2 0.98). A similar consideration for the effect of the 14 micron ozone band within the 14 to 16 

micron range is not as favorable. In the case of the standard atmosphere, for an h 0 of zero, the ozone 
diagram for 14 to 16 microns shows a transmissivity from the surface of only about 0.8. However, 
the emission originates mostly from the 10 to 25 km low temperature range, which also dominates the 
carbon dioxide diagram. The net effect of the residual 0 3 emission is then the same as the C0 2 emis- 
sion, and it was therefore neglected. 

The window from 10.75 to 11.75 is so transparent that small amounts of solid particles, or drop- 
lets, in the air affect transmissivity. For example, cirrus clouds, often invisible from the surface, 
or dust trapped near the tropopause seem to have more influence on the outgoing radiation in the win- 
dow than commonly thought (Reference 14). This issue has been neglected, since very little is known 
about the altitude distribution of aerosols and dust. The effect, if considered, would be a reduction in 
the radiance emitted from the earth reaching satellite altitudes and a possible increase of radiance in 
the clear atmosphere for h 0 values between 5 and 20 km. 


RESULTS OF THE INTEGRATION 

The results of the integration are summarized, for each wavelength interval, in Figures 10 
through 14. In the diagrams, the radiance N (w/m 2 -ster) within the spectral band and the corre- 
sponding blackbody temperature (°K) are plotted versus altitude h 0 above or below the geometrical 
horizon (h 0 = 0). Included at the left side of each figure is the radiance value for r Q = 0, which is 
the radiance emerging from the earth in the direction of the radius vector. 

A comparison of the radiance for r Q = 0 (center of disk) and h 0 = 0 (horizon) shows a small but 
definite limb-brightening effect for all atmospheres in the 14 to 16 micron band (Figure 10). Almost 
all calculated model atmospheres show a maximum in the vicinity of +20 km, and an almost identical 
behavior above 30 km. The arctic winter atmosphere has the same general appearance but deviates 
somewhat in the absolute levels. The reasons for this are the low temperatures of the upper atmos- 
phere and the associated low densities, which are consistent with the low temperatures in the tropo- 
sphere. Both effects can be observed in Figure 5. Clouds have a noticeable influence on the radiance 
in the tropics only— and there, just near the center of the disk. In the strong carbon dioxide band, 
the earth appears very uniform with the maximum radiance near the rim and with a generally warmer 
summer and cooler winter pole. 

The appearance of the earth in the atmospheric window (Figure 11) is quite different from that 
seen in the carbon dioxide region; the slope of the horizon is generally much steeper. Only in the 
arctic winter case, the geometrical horizon becomes visible. In practice, this will occur rather 
rarely since, for the most part, clouds will obscure the view, making the disk of the earth appear 
rather irregular. Enormous changes in the radiance levels from one point of the horizon to an adja- 
cent one must be expected; for example, deserts emitting about 10w/m 2 -ster may be next to high 
clouds emitting only 0.9 w/m 2 -ster. The meteorological situation has a great influence on this 
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Figure 10-Curves of radiance in the 15 micron C0 2 band vs. h 0 ( r ) for the five model 

atmospheres. Solid lines represent clear sky, and dashed lines represent high cloudiness 
conditions. 
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SUMMER 240 ‘ 

ARCTIC WINTER-^: _ 

220 

ARDC 210 ~ v " T • 

STANDARD 200 

TROPICAL 


Figure 1 1 -Curves of radiance in the atmospheric window vs. h 0 (r 0 ) for the five model at- 
mospheres. Solid lines represent clear sky, and dashed lines represent cloudiness conditions 
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Figure 12— Curves of radiance in the 6.3 micron H 2 0 band vs. h fl (r„) for four model at- 
mospheres. Solid lines represent clear sky, and dashed lines represent high cloudiness 
conditions. 



Figure 13-Curves of radiance in the rotational H 2 0 band vs. h 0 (r 0 ) for four model atmos- 
pheres. Solid lines represent clear sky, and dashed lines represent high cloudiness conditions. 



Figure 14— Curves of radiance in the 9.6 micron 0 3 band vs. h 0 (r 0 ) for the tropical Arosa 
normal, Arosa high, and arctic winter atmospheres. Solid lines represent clear sky, and dashed 
lines represent high cloudiness conditions. 


wavelength region. The threshold of a horizon detector system must be set low enough that high 
cloudiness will not be mistaken for outer space. 

The 6.3 micron water vapor band shows less contrast than the window at the center of the disk, 
r 0 = 0, (Reference 15) and even less near the horizon (Figure 12). The low temperatures at the 
tropical tropopause and an increase in the water vapor mixing ratio in the stratosphere give rise to 
a warm zone at 25 km above a partial minimum near 15 km. The same effect creates a small plateau 
in the desert atmosphere between 20 and 30 km. Tropical clouds, however, are still visible near the 
horizon. 

The general shape of the horizon in the rotational water vapor band, shown in Figure 13, appears 
much more favorable, although the absorption coefficients are of the same order of magnitude as in 
the 6.3 micron band. The reason lies in the energy distribution of the Planck function. The radiance 
in the 6.3 micron band is proportional almost to the 10th power of temperature within the tempera- 
ture limits of 200° and 250°K, whereas the rotational water vapor band in this same region is propor- 
tional to a power of only 2 or 3. The contrast between radiance values across the disk of the earth 
is therefore much less than in the 6.3 micron band. In the water vapor band, all atmospheres show a 
steep decline at about a 10 km altitude. Tropical clouds again are visible, and minor humps occur 
similar to those in the other water vapor band at 6.3 microns. 

In the ozone absorption region, shown in Figure 14, the clear atmosphere exhibits a very pro- 
nounced limb-darkening effect. Similar to the window region, strong contrasts between clear and 
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cloudy areas exist. The stratification of ozone yields a flat minimum or a small plateau between 10 
and 30 km, which is analogous to the situation for tropical water vapor. We can even speak ol a double 
horizon, one near 40 km and the second one near 5 km— the latter one appearing only under clear sky 
conditions. A comparison of the Arosa ozone distributions shows a small but interesting phenomenon: 
An increase of about 50 percent in the ozone content in the upper atmosphere decreases the radiance 
of the disk (only by a small amount, perhaps 8 percent) but causes a rise in radiance for h 0 values 
above 5 km. The arctic winter atmosphere is also exceptionally low in radiance. Again, in the ozone 
region, clouds contrast strongly against clear skies. 


COMPARISON OF CALCULATED VALUES TO TIROS RADIATION DATA 

The Tiros meteorological satellite radiation experiment provides observational data for com- 
parison with some of the above calculations. One channel of the radiometer was sensitive to radia- 
tion in the 8 to 13 micron region, and it is therefore possible to make a rough comparison with this 
observed region by combining the calculations for the 10.75 to 11.75 micron and 8.9 to 10.1 micron 
bands. While the 5-degree field of view of the Tiros radiometer does not permit sufficiently high 
resolution of the horizon, it is possible to compare data taken with the sensor axis nadir angle near 
zero with the calculations made for r 0 = 0. 

It is assumed that the radiance per wavelength increment seen in the 8 to 13 micron region would 
be the same as that seen in the 10.75 to 11.75 micron region, with the exception of the 8.9 to 10.1 
micron ozone band, where the calculated radiance for that region was used. In this way the total 
blackbody temperature was derived. The temperature thus calculated for the standard atmosphere 
and clear skies was 280° K compared, for example, with an average of 289°K for clear areas over the 
eastern United States (Reference 16). The calculated temperature for the tropical atmosphere and 
clouds near the tropopause of 201 K compares favorably with the observed temperature over Hurricane 
Anna of about 204 K (Reference 17). 

The above comparisons show reasonable agreement, considering the strong dependence of the 
radiance on the temperature profile and gas distributions. A more thorough comparison of observa- 
tion and theory would necessitate soundings of temperature and gas distributions at the times the 
satellite data were acquired. Calculations could then be made, using the sounding data, and could be 
compared with the radiation observations. 


CONCLUSIONS 

For the purpose of horizon sensing, the choice of a spectral region is governed by a number of 
considerations. One important design parameter is the lowest value of radiance N min within the disk 
of the earth. A horizon sensor must identify this level as the earth, in contrast to space. Since in 
certain design concepts the detector receives radiation only from a small region near the horizon and 
never from the central part of the disk, the minimum radiance for h 0 = 0 is important. Table 3 lists 
a value N 0 slightly below the minimum radiance, which could serve as a conservative estimate of this 
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Table 3 


Summary of Important Parameters. 


Wavelength 

(microns) 

N = 0) 

(w/m 2 -ster) 

N (h 0 - 0) 

(w/m 2 -s ter) 

N „ 

(w/m 2 -ster) 

A h 
(km) 

N 

max 

N 

ma x 

N 

m i n 

N 

m i n 





fro = 0) 

(ho = 0) 

6.33 - 6.85 

0J 02 

0.201 

0.08 

29.2 

5.69 

2.05 

8.9 - 10.1 

1.09 

1.57 

0.9 

14.7 

8.95 

2.79 

10.75 - 11.75 

0.86 

0.86 

0.7 

13.4 

14.3 

10.9 

14 - 16 

3.54 

3.69 

3.0 

13.3 

1.61 

1.68 

21 - 125 

13.3 

14.7 

12 . J 

7.8 

1.74 

1.37 


threshold value. The difference Ah between the highest and lowest horizon, quoted for the threshold 
level n 0 , is important only for sensors capable of resolving the true shape of the horizon. 

The curves in Figures 10 through 14 may also be interpreted as time functions obtained from an 
ideal detector of zero time constant, without noise, and with an infinitely narrow field of view scanning 
across the horizon. A more realistic time- constant- to- scan speed ratio and a finite field of view can 
be considered to be bandwidth restrictions (electrical filters) applied to the Fourier spectrum of the 
horizon functions. The influence of detector noise on horizon determinations was discussed recently 
by Goetze (Reference 18). 

For many instrument designers the ratio between the maximum and the minimum value of N be- 
comes very significant; these ratios are listed for the center of the disk (r 0 = 0) and for the horizon 
(h„ = 0 ). 

The great advantages of the long wavelength, high absorption regions between 14 and 16 microns 
and beyond 21 microns are apparent. The radiance values are higher; uncertainties in the horizon 
Ah caused by meteorological and seasonal effects are lower; and the contrast across the disk of the 
earth is low in comparison with the values in the shorter wavelength regions including the window. 
Unfortunately, instrumental techniques for wavelengths longer than 30 or 40 microns are not as highly 
developed as those at shorter wavelengths. However, the 15 micron carbon dioxide band is well within 
the convenient operating range of thermal detectors, and filters are commercially available. 

The comparison of spectral regions in Table 3 includes seasonal effects as well as local mete- 
orological phenomena. It can be expected that the two effects enter accuracy considerations in a 
different manner. For example, the low arctic winter radiance in the 15 micron carbon dioxide band 
is semi- stationary and could, at least theoretically, be compensated for in the same manner as we 
could compensate for the equatorial bulge of the earth. Meteorological effects, on the other hand, are 
more or less random and cannot be compensated for, other than by integration over time and space. 
From these considerations, the 15 micron band becomes even more attractive. 
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Considering both the theoretical and instrumentational aspects of the problem, it is concluded 
that a scanner sensitive in the interval of approximately 14 to 16 microns is most promising for 
minimal errors in the determination of the local vertical from a satellite by means of a horizon sens- 
ing system. 

It is planned to include a 14 to 16 micron scanner on one of the Tiros satellites in the near future 
in order to compare theoretical with observed results in this spectral interval. 
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Appendix A 

List of Symbols 

Spectral radiant emittance (Planck function) for a blackbody of temperature T (watts/m’ - 
micron) 

Acceleration due to gravity 

Height above the surface of the earth 

Height of line s (along radius vector normal to s) 

Generalized absorption coefficient 
Infrared radiance (watts/m 2 -ster) 

Conditions of normal temperature T 0 and pressure P 0 
Atmospheric pressure 

Normal pressure (760 mm Hg, 1013.25 mb) 

Pressure of 10 mm Hg, above which the pressure correction factor P/P 0 for calculating the 
u* for ozone remains constant atP^P,, = 0.132 
Distribution function of an optically active gas 
Radius of assumed spherical earth (6370 km) 

Gas constant for dry air (2.8704 x 10 6 cm 2 /sec 2 - °K) 

Geocentric distance 

Geocentric distance to line s (length of radius vector normal to s) 

Straight line path along which the integration of u* is carried out 
Kinetic temperature 

Equivalent blackbody temperature corresponding to integrated radiance within a given spectral 
interval 

Normal temperature (0°C, 273.16°K) 

Actual optical thickness of an absorbing (and radiating) gas 
Reduced optical thickness, corrected for pressure and temperature 
Zenith angle of s 
Wavelength (microns) 

Density of air 
Density of air at NTP 

Transmissivity for straight line beam radiation 
Spectral response of a sensor 
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